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Abstract 



In this paper, we establish rigorous existence theorems for a mathematical model of a thin 
inflated wrinkled membrane that is subjected to a shape dependent hydrostatic pressure load. We 
are motivated by the problem of determining the equilibrium shape of a strained high altitude 
large scientific balloon. This problem has a number of unique features. The balloon is very thin 
(20-30 yum), especially when compared with its diameter (over 100 meters). Unlike a standard 
membrane, the balloon is unable to support compressive stresses and will wrinkle or form folds 
of excess material. Our approach can be adapted to a wide variety of inflatable membranes, but 
we will focus on two types of high altitude balloons, a zero-pressure natural shape balloon and a 
super-pressure pumpkin shaped balloon. We outline the shape finding process for these two classes 
of balloon designs, formulate the problem of a strained balloon in an appropriate Sobolev space 
setting, establish rigorous existence theorems using direct methods in the calculus of variations, 
and present numerical studies to complement our theoretical results. 

1 Introduction 



Balloons play an important role in NASA's current scientific investigations, including upper atmo- 
sphere research, high energy astrophysics, stratospheric composition, meteorology, and astronomy. 
With the development of the Ultra Long Duration Balloon (see, e.g., and the possible uses of 
balloons in the exploration of planets in our solar system, balloons will play an important role in 
NASA's future scientific endeavors (see, e.g., j22j)- Furthermore, many of the techniques that have 
been developed for the analysis of balloons can be readily adapted to other light-weight membrane 
structures, including solar sails, inflatable rovers, aerobots, and gossamer spacecraft. 

Large scientific balloons are regularly flown by NASA to carry out research in the stratosphere. 
With a fully inflated diameter of over 100 meters and a thickness of 20-30 microns, a large scientific 



balloon does not behave like a standard inflated membrane. Because the balloon is so thin, it is 
unable to support compressive stresses and will instead form folds of excess material or wrinkle. 
Although our theoretical results could be applied to partially inflated balloons with the contact 
problem handled by constraints, for the numerical simulations considered in this paper, the balloons 
will be fully or near fully inflated, so wrinkling can occur, but folds cannot form. Because the diameter 
of a typical balloon is so large in comparison with the thickness of the balloon skin, the problem 
of a strained large scientific balloon falls outside the realm of typical inflated membranes. For a 
survey on inflated membranes, including wrinkling, see |23l Chapter V, Section T] and the references 
therein. To handle wrinkling, one replaces the wrinkled region by a smoothed out pseudo-surface. 
This can be done via pseudo-constitutive relations as in the work of Stein and Hedgepath |27j where 
a variable Poisson ratio is introduced. We follow the approach of Pipkin [23] where a "relaxed" strain 
energy is used to model wrinkling (see also |26j). In addition to wrinkling, the balloon problem is 
characterized by a number of other features, including large displacements, relatively small strains, 
and a shape dependent hydrostatic pressure load. Problems in nonlinear elasticity have been well- 
studied using a variety of methods (see, e.g., [2], [23] ), including direct methods in the calculus of 
variations (see, e.g., ^H], [IB Appendix A]), but our existence results for a thin inflated strained 
large scientific balloon are the first rigorous analytical treatment of this class of membrane problems. 
Our analytical results affirm the use of numerical models based on optimization, and help explain 
the efficacy of these types of schemes (see, [S]-[E])- 

To date, the workhorse of NASA's large scientific balloon program has been the zero-pressure 
natural shape balloon, a design that goes back to the 1950's (see |T]). In recent years, due to 
scientists' demands for long duration mid-latitude balloon flights, a design concept that has come to 
be known as the pumpkin balloon has been in development. While the design shape of a balloon is 
a theoretical target which the fully inflated balloon is intended to assume, shape generating models 
are usually limited to weight and pressure considerations and ignore straining in the film. However, 
upon inflation, the real balloon envelope responds to the differential pressure and strains. For these 
reasons, the analysis of an inflated balloon usually is divided into two distinct parts, shape finding 
and stress analysis. Shape finding is carried out by the balloon designer. Typically, the fully inflated 
balloon must maintain a payload at a constant altitude. At float conditions, the shape of the 
balloon and its volume must be determined. Archimedes' principle states that the upward force 
(lift) generated by the lifting gas is equal to the total weight of the balloon system. At equilibrium 
the lift is equal to the weight of the air displaced by the lifting gas less the weight of the lifting 



gas, i.e., Lift = gp a ^o — gPg^o, where g is the gravititational constant, p a is the density of the 
atmosphere, p g is the density of the lifting gas, and ujq is the volume of the gas bubble. p a and p g 
are functions of altitude, temperature, etc. and are assumed to be known. The specific buoyancy is 
defined to be b = g(p a — Pg)- To emphasize that a quantity depends on the design conditions, we 
add a subscript of l d\ For example, at float, we have Lift = b^ujQ^. We will assume that for other 
altitudes corresponding to specific buoyancy b and volume ojq, the following relation holds 



Note, b is decreasing as a function of altitude, so means that volume is increasing as a function 
of altitude. For example, a balloon with a sea level volume of ojq will expand to approximately 300wo 
at 39 kilometers. 

The shape finding processes for the two designs discussed here are outlined in Section 2. More 
detail on shape finding can be found in 4 . A byproduct of the shape finding process is a gore 
pattern that the manufacturer uses to fabricate the balloon. A gore is a long flat tapered panel of 
film. The gores are sealed edge-to-edge in such a way that when fully inflated, the complete balloon 
will assume a shape very similar to the one desired by the balloon designer. The shape of the gore 
pattern is an important input into the stress analysis. 

The second part of the balloon problem involves an analysis of a pressurized elastic membrane. For 
efficiency, the balloon needs to be as light as possible, yet it must be strong enough to operate safely 
over its service life. A meaningful stress/strain analysis requires knowledge of certain mechanical 
properties of the structural elements (e.g., film thickness i, Young's modulus E, and Poisson's ratio 
v). The membrane is assumed to be made of a linearly elastic isotropic material, and because it 
is so thin, it is unable to support compressive stresses, and will wrinkle instead. Following the 
approach of Pipkin [21], we model wrinkling by relaxation of the film strain energy density. Pipkin's 
approach was adapted to large scientific balloons by Collier (see, e.g., ^Sj) and has been used in 
a number of subsequent papers (see, e.g., jSJ-El)- We assume that the balloon is in a quasi-static 
steady equilibrium state and because all of the forces on the balloon are conservative, the problem 
of determining an equilibrium shape can be formulated in terms of a variational principle. This 
formulation lends itself to theoretical analysis and efficient numerical computation. While the main 
results of this paper are the existence theorems in Section |IJ we also include numerical solutions. 

This paper is divided into six sections. In Sectional we provide some historical background on 



buj = b d u 0jd . 




mathematical models that are used for the shape finding process. The most common models are 
derived from the equations for an axisymmetric membrane. We also consider a cyclically symmetric 
pumpkin balloon. In Section |31 we present a mathematical model for a strained balloon shape. 
The membrane is modeled as a nonlinearly elastic membrane shell (see, |121 Sec. 9.4]). It will be 
important to distinguish at least three types of balloon configurations: (a) Q C IR 2 - the flat reference 
configuration or natural unstrained state; (b) Sq C IR 3 - an initial curved configuration; (c) S C IR 3 - 
a general configuration (not necessarily in equilibrium) . f2 and So are generated by the shape finding 
process and S is generated in the analysis of a strained shape. Our numerical solution process is 
robust and So need not be close to the final strained equilibrium state. 

Many non-shallow shell theories start with an initial curved reference configuration in equilibrium 
for some load, and formulate the equilibrium equations in terms of an appropriate displacement field 
u from that state corresponding to a load increment. Since the actual balloon is constructed from a 
number of thin flat tapered panels of film, it is more natural to parameterize our problem in terms of 
displacement from the flat reference configuration, i.e., x : — > IR 3 . This is particularly important 
when considering pumpkin balloons, where any curved configuration is necessarily strained and/or 
wrinkled and there is no obvious equilibrium state So to begin the solution process. In our approach, 
we start with a shape So (usually obtained from the shape finding process), but do not require it to 
be in equilibrium. Our film strain energy density function Wf is equivalent to the two-dimensional 
strain energy of Koiter's nonlinearly elastic membrane shell (see, |12l p. 450]). We will assume the 
balloon is a linearly elastic isotropic material. Our finite element model (FEM) triangulates f2, then 
uses piecewise linear elements and constant strain triangles in its implementation. As we will see, 
closed subspaces of W ' {Q) are the natural setting in which to study the balloon problem. In order 
for a mathematical model to produce meaningful results for this type of membrane, it is important to 
take wrinkling into account. This can be accomplished by replacing the usual strain energy density 
Wf with Wf, the quasiconvexification of Wf (see |24|). Wf is the largest convex function that does 
not exceed Wf, and in our case, we are able to calculate an explicit formula for Wf. Using the 
relaxed strain energy density of the membrane, the existence of equilibrium balloon shapes follows 
from direct methods in the calculus of variations (see, e.g., [16 J. 

The problem of determining the strained state of a balloon is formulated in Section|3]and existence 
results are established in Section^ Variations on the model presented in Section 3 have been used 
for numerical studies of strained balloons (see, [3]- [5]), however, Section 0] contains the first rigorous 



existence results for a model of this type. The tendon model implemented in earlier work (see, e.g., 
EUGi]) treated the tendons as linearly elastic strings whose strain energy was added directly into 
the total energy of the balloon system. In the present paper, we focus on inextensible tendons, so 
tendons are included through the use of constraints. Using a very stiff tendon in the models 
we obtain results that are consistent with the results presented here. To complement our theoretical 
results, we will present numerical solutions of strained pumpkin and strained zero-pressure natural 
shape balloons in Section |SJ Our theoretical results apply to asymmetric balloon shapes as well 
as cyclically symmetric ones, but we will limit our numerical studies here to cyclically symmetric 
shapes. See [7] for examples of asymmetric balloon shapes. Section includes concluding remarks. 

2 The shape finding process 

Before moving to the main results of this paper, it is important to have a clear picture of the shape 
finding process and how it determines a gore pattern. We focus on the axisymmetric zero-pressure 
natural shape balloon and the cyclically symmetric pumpkin balloon. 

2.1 Zero-pressure natural shape balloons 

The design shape of a high altitude balloon is normally based on conditions at maximum (i.e., float) 
altitude. In addition to known quantities such as film weight density per unit area uu, suspended 
payload L, and specific buoyancy b, there are other parameters relevant to the design shape, includ- 
ing the circumferential stress a c and a constant differential pressure term po (po is the differential 
pressure at the base of the balloon, Z = 0). The usual force balance analysis for an axisymmetric 
membrane involving weight and differential pressure (see, e.g., [21 p. 343]) leads to a system of 
ordinary differential equations 



for the generating curve T(s) = R(s)i + Z(s)k, where R'(s) = sin#, Z'(s) = cos 6, a m is the 
meridional stress, f = —ph — iok, p = bZ + po is the differential pressure, t is a unit tangent to 
the generating curve, b is the inward normal, and s is arc length of the generating curve. See 
Figure HJa). t makes an angle of 6 with k. The length of the generating curve and 6(0) = 9q 
are unknown at the start, and are normally found by using a shooting method to solve a nonlinear 
boundary value problem involving (|'2.1|) and auxiliary conditions (e.g., the total meridional tension 
in the k direction at the bottom of the balloon must equal the weight of the suspended payload, 
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Figure 1: T(s) = R(s)i + Z{s)\l; (a) Zero-pressure natural shape gore generator; (b) Pumpkin gore generator T. 
Unit vectors t, b, i, and k are indicated. 

Lj cos 6*o = 27r(i?c" m )(0); zero weight at the top of the balloon implies R(Ld) = 0). See [I] for further 
details. 

The model that is used most commonly for applications to large scientific balloons is the zero 
pressure natural shape (ZPNS) model which assumes po = 0, a c = and all tension is carried in the 
meridional direction. While the assumption of zero hoop stress is clearly violated in a strained float 
shape, the ZPNS design has proven to be successful for NASA. The original work on high altitude 
plastic balloons was carried out in the 1950's at the University of Minnesota, where the term natural 
shape balloon emerged (see, [Q). 

The shape finding process for the ZPNS balloon leads to an approximation of the axisymmetric 
shape by a cyclic shape where a fundamental region Qp C 1R 3 is taken to be a developable surface 
(see Figure Ufa)). Thus, there is an isometry that takes Qf into a plane, from which the "lay- 
flat" gore pattern Gf is determined. Gf is used by the manufacturer in the construction of the 
balloon. A typical large zero-pressure natural shape balloon has over 100 gores and it would be hard 
to distinguish the axisymmetric shape from the cyclically symmetric shape generated by the Qf- 
However, in Figure [2fb), we use a small number of gores in order to highlight certain features. 

2.2 Pumpkin balloons 

In the past, long duration balloon flights took place in regions such as Antarctica where the wind 
currents are regular and the balloon can circumnavigate the South Pole under the same thermal 



conditions for a long period. In 2005, for example, the Cosmic Ray Energetic and Mass Experiment 
(CREAM) using a ZPNS balloon flew for nearly 42 days, maintaining an altitude between 38 and 
39 kilometers for most of its flight. However, at mid-latitudes, a significant amount of ballast must 
be carried by a ZPNS balloon in order to maintain constant altitude over several diurnal cycles. 
At night when the temperature of the lifting gas cools, ballast must be dropped, and during the 
day when the temperature of the gas rises, gas must be vented. Typically, a balloon carries enough 
ballast for a few of these cycles, restricting the length of a ZPNS mid-latitude flight to no more than 
a few days. One way to avoid the need to carry significant ballast is to construct a balloon that 
can contain a sufficient amount of gas to maintain altitude at night and is strong enough to hold 
the overpressure caused by solar heating during the day. This led to the consideration of a balloon 
design that has come to be known as the pumpkin balloon, a term coined by Smalley in the early 
1970's (see, |25j). The pumpkin shape was fore-shadowed by Taylor in his work on parachutes (see 
|28l Fig. 1]). Currently, NASA's Ultra Long Duration Balloon (ULDB) Program aims to develop 
a pumpkin balloon that is capable of staying aloft for one hundred days at any latitude. As of the 
writing of this paper, the ULDB pumpkin has exhibited deployment problems, but these issues are 
outside the scope of this paper. For more on this topic, see 0-0 • 

The principal behind the pumpkin balloon is to use a light-weight film as a gas barrier and strong 
reinforcing tendons for pressure confinement and to carry the weight of the balloon system. Roughly 
speaking, increasing the curvature in the circumferential direction has the effect of transferring most 
of the load to the tendons. To model the lobing in an ideal pumpkin gore, one can assume a 
fundamental region Qp is a subset of a tubular surface. Consequently, the shape Qp is doubly curved 
(see Figure |2| . If a pumpkin gore is constructed from a flat panel of film, any inflated pumpkin 
configuration is necessarily strained and/or wrinkled. 

Next, we define pumpkin balloons that are analyzed in this paper. We begin with a description 
of a constant bulge radius pumpkin that is parameterized as a tubular surface. Let T(s) = R(s)i + 
Z(s)k G IR 3 be a planar curve that we call the generator of the pumpkin gore. Arc length measured 
along the generating curve is denoted by s. See Figure ^b) for a sketch of a representative T drawn 
in the RZ -plane. Apriori T is unknown, and must be derived from equilibrium conditions. A detailed 
exposition of the shape finding equations for a pumpkin balloon, including the determination of T, 
is presented in [4 . In the following, we will assume that T is known. The generator is parameterized 
by s and a prime indicates differentiation with respect to s. Let t denote the unit tangent and 
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(a) ZPNS Q F (b) Pumpkin Q F 



Figure 2: (a) Qf for a zero-pressure natural shape balloon; (b) Tubular surface Qf with generating curve T and 
bulge radius tb- 

b = t x j the inward unit normal of T; 9 = 9{s) is the angle between t and k (see Figure ^b)). The 
set {b,t,j} gives a right-hand curvilinear basis for K 3 . The curvature of T is denoted by k where 
T"(s) = «(s)t'(s). We define a tubular surface with generator T and constant radius rs as follows 
(see, [HI p. 89]): 

x(s, v ) = T(s) + tb (j shif — b(s) cos v) , \v\ < vb(s), < s < Lj, (2.2) 

where vb(s) is known from the shape finding process and depends on the number of gores N. See 
Figure ^for a typical generating curve and Figure Efb) for Qp as determined by the shape-finding 
process for a pumpkin balloon based on a tubular surface model. 

We assume that the pumpkin gore Qp is situated symmetrically with respect to the y = plane 
and interior to the wedge defined by the half-planes y = ±t&n(ir/N)x where N is the number of 
gores and x > 0. We will refer to tb as the bulge radius of the pumpkin gore. The curve traced by 
v — > T(s) + 7\b(j sin w — b(s) cos v) is a circle lying in the plane with normal t(s). By construction, 
the length of the circular arc that is formed in the pumpkin gore is 2rBVB(s) (see Figure Hfb)). We 
call vb(s) the bulge angle. For large balloons with a large number of gores, < tbk{s) « 1, but in 
Figure Hfb), we used a small number of gores in order to highlight certain features. A typical large 
pumpkin balloon has over 200 gores, yet the lobing in a fully deployed pumpkin is clear to see (see 
[HI Figure 2(a)]). We define the theoretical three-dimensional pumpkin gore Qf to be 

Qf = {(x,y,z) = x(s,v), -v B (s) <v < v B (s), < s < L d } . 



(a) Qf (b) Gf - lay-flat of gore 



Figure 3: (a) £/F-initial three dimensional gore configuration; (b) Qf = Gf - usual lay flat component for one 
pumpkin gore. 

A complete shape So has cyclic symmetry and is made of N copies of Qf- Note that the length of 
the centerline of Qp is 



respective centerlines of Qp and Gp are isometric. The length of a rib (i.e., a circular arc) in the 
spine of Qf is 2r^us(s) and this is the same as the length of a corresponding segment orthogonal 
to the centerline of Gf. It follows that L s , the corresponding edge length of the nominal lay- flat 
pattern Gp, is longer than the tendon length Lj. For the pumpkin example considered in Section 5, 
L s = 1.005L£. A tendon is normally encased in a sleeve that is sealed along the length of a gore seam. 
To accommodate the lack-of-fit, the sleeve is gathered first. The tendon is tacked to the gathered 
sleeve at a number of locations. The amount of gathering is not uniform. Maximum gathering takes 
place near the equator, and decreases as one gets closer to the central axis of the balloon. The 
smoothed out sleeve and the edges of two adjacent gores are heat sealed together to form one seam. 




Normally, the length of a tendon is taken to be 




Corresponding to Qf C IR 3 



is the lay-flat configuration Gf C M 2 shown in Figure [5(b). The 
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Figure 4: (a) Q.f; seam is parametrized by (h(v),v) for < v < L c ; (b) close up of gore bottom; IV component of 
Qf that will be attached to the end fitting. 

The interested reader is referred to jU] and the references therein for a more detailed discussion of 
tendon-foreshortening and how it affects the film stresses. 

3 A model for a strained balloon 

In this section, we formulate the problem of a strained balloon and describe the components in our 
variational principle and related constraints. 

3.1 Preliminaries 

We first define the reference configuration for a complete balloon. We begin with a discussion of 
the lay-flat configuration that is obtained from the shape finding process. Let £If = Gf and let L c 
denote the gore length where 



and h(v) and L c are known from the shape-finding process. We can assume that u = h(v) is a 
function, so that is a C 0,1 domain. Here, Gf is either a ZPNS gore or a pumpkin gore, but in 
other applications it may be some other reference configuration. In terms of h(v), the seam length 
of fip is (see Figure Efb)) 



Q F = {(u,v) 



h(v) < u < h(v), < v < L c } 




Let 



I> = {(u, v) | - h(0) <u< h(0),v = 0} . 
Y*F = {(u,v) I u = h(v), < v < L c } . 

Let A > 0, f2i = &f, Ti = Tp, Si = S^, and h\(v) = We define hi{v) = /ij_i(t>) + A, 

^ = {(n + (i - 1)A, f ) | (u, v) G ill} 
T 4 = {(u + (i-l)A,v) | (u,«)€r 1 } J 
S, = {(u + (i-l)A,v) | («,w) G Si}, 

where i = 2, . . . , N. The balloon is constructed from N identical flat gores Oj that are sealed edge 
to edge (in practice, they overlap slightly and 150 < N < 290). The right seam of is identified 
with the left seam of fit+i. The right seam of O^v is identified with the left seam of Q±. With these 
identifications, we define 

N 

n = |J n f . (3.i) 

Remark 1 Due to the slight overlap of the edges of the fij's as described above, we can think of 
£1 as defining a two dimensional manifold. For convenience of exposition, we can extend slightly 
the gore width of the fij's (call the extended gore so that adjacent local charts overlap (i.e., 
Cli n Sj n fii+i 7^ 0) and has a manifold structure. Since the overlap is so small in our applications, 
we treat it as a curve S,. The metric is the usual Euclidean metric on each fli C IR 2 . Open sets 
can be defined using this metric with special attention given to points that lie along the seams, 
i.e., (u, v) G Sj C fi. In what follows, we study the balloon problem in a certain Sobolev space. 
While it is possible to extend the necessary definitions and results such as the Sobolev Embedding 
Theorems to a manifold setting (see, e.g., [3]), introducing additional notation for this purpose would 
be cumbersome. Since it only makes sense to consider the differentiability of x : f2 — ► IR 3 on the fij's, 
we adhere to the following conventions. We write x £ C 1 (r2,IR 3 ) to mean x £ C l (Q \ U^iS^IR 3 ). 
With an appropriate metric defined on CI, we can consider x G C(il,IR 3 ), where the boundary of Vl 
is T = \jf = iTi. With these conventions, we consider mappings x G C 1 (fi,IR 3 ) n C(Cl, K 3 ). 

Typically, there is an end-plate located at the top and bottom of a balloon, but for simplicity, 
we assume that a gore comes to a point at the top. T corresponds to the bottom end-plate. For ease 
of exposition, we will assume that the material properties (e.g., Young's modulus, Poisson's ratio, 
film thickness, etc.) are constant over Cl. We will formulate the equilibrium equations for a complete 



balloon, but in practice, we will often assume symmetry in order to reduce the number of degrees 
of freedom in our numerical model. We assume that dp = Q\ is situated in the w-plane with its 
center along the f-axis (see Figure 2}. 

A deformation mapping for a complete balloon is given by 

x : n -» S C IR 3 , 

where S = x(Tl) represents a balloon configuration and x G T> where T> is the set of admissible 
deformation mappings. Later, we will give a precise definition of T>. Let Sq = xq(£1) denote the 
initial configuration of the deformed balloon in 1R 3 . So need not be in equilibrium. In our formulation 
of the problem, the boundary of S has one component (i.e., xq(T), the part that is attached to the 
end-plate). For any x 6 T>, we require 

x(u, v) = Xq(u,v), (u,v) £ T. 

For future reference, we note 

cti(v) = x(hi(v),v) for < v < L c 

parameterizes x\^ t C IR 3 , the deformed seam between and Oj+i for i = 1, ...,N — 1. ctjv 
parameterizes the deformed seam between f2jv and fii. 

We define 

\ x \i,p = Mi,p + \v\i,p + l z li,p' 

where x(u, v) = x(u, v) ■ i, y(u, v) = x(u, v) • j, z(u, v) = x(u, v) ■ k, 

= / \ W \ P dA + [ \Vw\ p dA, 

\X7x\ 2 = tr (\7x T \/x s j = \x u (u, v)\ 2 + \x v (u, v)\ 2 , and \x\ = y/x • x. '•' denotes the usual Euclidean 
inner product in IR 3 . We will follow the convention that x \ = x u and x 2 = x v . Let 

W 1 ' p (n,M 3 ) = {x : |x|i,p < 00}. 



Next, we define an admissible deformation mapping. 



Definition 1 Let xq{u,v) = xq(u, v)i + yo(u, v)j + zo(u,v)k E W 1,p (fl, IR 3 ). An admissible defor- 
mation mapping, x(u,v) = x(u,v)i + y(u,v)j + z(u, w)k is a mapping such that 

(a) x E ^(fljn^fi); 

(b) x(u, v) = xq(u, v), (u, v) E T; 

(c) x(ui,L c ) = y(ui,L c ) = 0, for Ui = (i - 1)A, i = l,2,...,N; 

z(ui,L c ) = z(uj,L c ) for all 

The set of all admissible deformation functions is denoted V. 
In the following, we let 

X = {x = x + x | i G Wo' p (0,IR 3 )} . (3.2) 

X is the completion of T> with respect to the norm |£c|i, p . In our formulation of the balloon problem, 
p = 4 and a solution is x = x + xq E X where a: E VFq' (f2,IR 3 ), and Wq IR 3 ) is the completion 
of C 1 (r2) fl Co(O) with respect to |sc|i,4. We will consider closed subspaces of X in the form X s , 
where X s is the completion of 

V s = {x E V | g(it,v,a;,Va!) < 0} 

with respect to |£c|i,4- g = (<7i,<?2, ■ ■ ■ >5m) will incorporate boundary conditions and constraints of a 
certain type that will ensure X s is closed. 

In the following, we let x\ B denote the restriction of x E T> to B C f2. Typically, £> will denote a 
curve in Q corresponding to the middle of a gore or a gore seam. Since p = 4 in our application, the 
trace of x E X for such a set i3 is well defined (see p. 257]). In particular, let x k ,x E X and 
\xk — £c| i,4 — > and consider g^s of the following types: 

1. Global Constraint 

gi(x)=0, for g t E C(X,TR). (3.3) 
From (|3.3|) . we see lim gi{x k ) = gi{ lim Xf~) = g%{x). 

k— >oo k^oo 

2. ZocaZ Constraint 

9i (x\ B ) < 0, for ^(^Ib) E C(0,IR). (3.4) 
From (J33J), we see lim gi{x k \ B ) = lim x k \ B ) = gi{x\ B ). 

k— >oo fc^oo 



We are led to the following. 

Lemma 1 Let X = [x = x + x I x e W 1,4 (fi, iR 3 )} and g = (gi, . . . ,g m ). If gi satisfies kH.ty) or 
for B = Bi C Q, i = 1, 2, . . . , m, then X s is a closed subspace, where 

X g = {a: = x + x | x £ Wq A (Q, ffi 3 ), g(u, v,x,Vx) < j . 

The proof of Lemma ^ follows directly from the definition of X and (|3.3[) - (|3.4j) . 

We will demonstrate that the components of g incorporating boundary conditions, symmetry 
conditions, or tendon constraints satisfy (|3.4j) and the volume constraint satisfies 1)3.3(1 . We first 
consider boundary/symmetry conditions. Tendon constraints are considered in Section f3, 3 . 1 1 and the 
volume constraint is considered in Section [3.3.21 

In order to keep the number of degrees of freedom in our numerical model manageable, we use X s 
in our simulations. Often, we consider a cyclically symmetric balloon, in which case, we only need to 
model one-half a gore, i.e., y = is a plane of reflectional symmetry and the right deformed seam of 
0,i must lie in the half plane y = tan(ir/N)x with x > 0. Let us suppose that g includes boundary 
conditions for a cyclically symmetric shape with N gores. In this case, there are two constraint 
boundary conditions. Let B\ = {(u, v) \ u = 0, < v < L c } and £>2 = Si, 

gi(v,x(u,v)) = x(u,v)-j, (u,v)eBi, (3.5) 
g 2 (v,x(u,v)) = x(u,v) ■ (j - tan(vr/iV)i) , (u,v) £ B 2 ., (3.6) 

gi = is the condition that ensures x(Qi) is symmetric with respect to the y = plane and g 2 = is 
the condition that cc(Si) lies in the plane y = tan(7r/iV)x. Suppose x^^x G X, \x^ — x\\^ — > 0, and 
g(xk) = for all k where g is defined in ()3.5|) - ()3.6|) . To demonstrate that X s is closed, we need to 
show g(x) = 0, but this is a consequence of the fact that g\ and g 2 are local constraints in the form 
(|3.4I) . A similar argument will apply to a lobed shape with dihedral symmetry D q . One boundary 
condition is in the form (|3.5|) . If we consider a periodic shape with q lobes and p gores per lobe 
where p, q are positive integers and pq = N, the second boundary condition implies that x(T, q ) lies 
in the y = t&n(ir/q)x plane, i.e., g 2 = x(u,v) ■ (j — tan(7r/g)i) = for (u, v) S H g . In this case, no 
geometric constraints are placed on tendons interior to the fundamental region. 

Remark 2 Although the equilibrium state of the balloon envelope will not have self-intersections, 
we do not put this restriction on x(£l) for x £ P. However, we do require x(£l) for an equilibrium 



state to be free of self intersections. A common theoretical and numerical approach to solving a 
membrane problem is to begin with an equilibrium configuration defined by a mapping, say y . A 
small increment is then added to the load and the corresponding displacement v 1 is determined 
for the new equilibrium, y + v 1 . The process is continued, and at each stage, v k is computed 
for k = 1, . . . ,n. The final solution is y + Y^k=i v ■ Because each intermediate solution must be 
an equilibrium solution, certain kinematical constraints are placed on the v k, s. In particular, each 
displacement field v k cannot give rise to an equilibrium state with self intersections. If a solution 
is desired whose load is much larger than the one corresponding to y , the number of intermediate 
steps could be large. By allowing self-intersections in the class of admissible displacements, we are 
able to proceed in a more direct fashion to the desired solution. This has the advantage of reducing 
the computation time. Although we allow intermediate configurations with self-intersections, in the 
end, we require that the equilibrium configuration is free of self-intersections. Note, a final solution 
could have regions of self-contact (see, e.g., [7j), but this could be described via local constraints in 
the form (|3.4f) . In this paper, the balloons are fully deployed so the contact problem does not come 
into play. 

3.2 Mathematical formulation 

Although film and tendon weight considerations are relevant to the shape finding process, they are 
not significant factors in a stress analysis for the applications in this paper, and for this reason, the 
film and tendon weight gravitational potential energies are not included here. We define the total 
potential energy of a balloon configuration for a deformation mapping x by Et = £t{x,Vx), where 

£t = E p + S f , (3.7) 

Ep is the hydrostatic pressure potential due to the lifting gas and Sf is the balloon film strain energy. 
Next, we discuss the properties of Ep and Sf. 

3.2.1 Hydrostatic pressure potential 

We follow the convention for hydrostatic differential pressure P that —P{z) > means that the 
inside of the balloon is pushing outward at height z units above the base of the balloon. For a 
super-pressure balloon, — P(z) = bz + po where b > is the specific buoyancy and po > 0. Typically, 
Po >> bztop- In a ZPNS balloon, po = 0. In our choice of coordinates, we assume that the base of 



the balloon is fixed and corresponds to z = 0, so that the potential for hydrostatic pressure P(z) is 
(see PU) 

E P = f P(z)dV = - [ (bz + p )dV, (3.8) 

JD JD 

where D is the region occupied by the gas bubble. Using the divergence theorem, (|3.8|) can be 
replaced by 

E P = - j {\bz 2 + p z)k ■ n dS, (3.9) 

where n is the outward unit normal to S, dS is surface area measure on S, and dD = S. The 
pressure distribution in an open system (i.e., without volume constraint) can be expressed as 

P(x) = -(bx-k + po). (3.10) 

Note that while b and po are known, the pressure potential depends on D. We can generalize the 
discussion by allowing a pressure distribution in the form P(x) = div(A(aj)), but for ()3.1U|) . we see 

A(x) = - (lb(x ■ k) 2 + Po (x ■ k)) k. (3.11) 

Note, A is not unique. If P is given by ([3.10^ . then by the divergence theorem, we have 

E P (x,Vx) = [ P(x) dV 

JD 

= [ f P (x,Vx) dA, (3.12) 
Jn 

where 

f P (x,X7x) = g P (x,&dj 2 X7x), (3.13) 

g P (x,A) = -{\b{x- k) 2 +p (x- k))k-A. (3.14) 

By |16l p. 117, Theorem 1.5], it follows that fp(x, ■) : IR 3,2 — > IR is polyaffine. However, fp polyamne 
implies fp and — fp are polyconvex. Hence, fp is quasiconvex ^3 p. 97]. 

Next, we establish a few inequalities. Since 

|adj 2 Va; • k| = \x t i x x )2 • k| < | Vcc| 2 , (3.15) 

it follows that 

- (\b\x\ 2 + |po| |sb|) |Vs| 2 < g P {x,<%<l} 2 Vx) < (±b\x\ 2 + \p \ \x\) \Vx\ 2 . (3.16) 
If \x\ < R, then 

- (^bR 2 + \po\R) \Vx\ 2 < g P (x, adj 2 Vcc) < [\bR 2 + \po\R) \Vx\ 2 . (3.17) 



Remark 3 Note, for a fully inflated shape po > 0, but if the balloon is not fully inflated, it is 
possible for po < 0. 



3.2.2 Balloon film strain energy 

Using methods in asymptotic analysis, Ciarlet derives the two-dimensional Koiter equations for a 
nonlinearly elastic shell in [121 Chapter 10]. The total energy of a shell of thickness 2e for an 
appropriate displacement field u is 



u 



I E af3ar,e ( a<Tr ( u ) _ 0(jt ) (a a fl(u) - O a fl) y/ady 
3 

+ L. f E a ^ T ' £ (b aT (u)-b ar )(b af 3(u)-b a p)V^dy, (3.18) 



where A e and /x e are the Lame constants, and E a ^ aT,e is the two-dimensional elasticity tensor of an 
isotropic shell, 

A\ £ ii £ , \ 

A e + 2/i e V / 

When applied to the balloon problem, we can ignore the second term in (|3.18|) . because e 2 is extremely 
small, i.e., we can ignore the bending or flexural energy. In our problem formulation, the natural 
state is the flat reference configuration. Thus, a a p = 5% where 5\ = 5% = 1 and Sf = 5\ = 0. We can 
express the strain energy in terms of a deformation from the natural state, i.e., 



x 



£ J E a P ar,e ^ . ^ _ ^ ^ . ^ _ ^ ^ ( g 2Q) 



recognizing that 7 a /3 = | ^a; Q • xp — dp) is the Cauchy-Green strain. Note, we do not linearize the 
Cauchy-Green strain in terms of the displacement field. 

The Lame constants are related to the Young's modulus E and Poisson's ratio u, via 

Ev fj^X + 2 M ) 

(1 + ^1-2^)' \ + » ' (321) 

M = 2(1 + !/)' 17 = 2(A + M ) 5 

where we have dropped the dependence on e, since we assume that the physical constants have been 
determined for a shell of fixed thickness. Previous work on balloons (see, e.g., |5]-|H]) utilized the 
expression in (|3.20j) for the total strain energy of the balloon membrane. For the convenience of the 
reader, we derive an expression for (|3,20|) that is consistent with the formulation in (SJ-jHl- 

We write the film strain energy Sf in the form 

S f = J W f dA, (3.22) 



where 

W f = \n : 7. (3.23) 

In (|3,23|) . n represents the Second Piola-Kirchoff stress tensor, 7 represents the Cauchy-Green strain 
tensor, and ':' is the tensor inner product. The contravariant components of n are denoted by n a @, 
the covariant components of 7 are denoted by 7 a( g, and n : 7 = n^ha- Assuming a linear elastic 
isotropic material, we have 

n a/3 = E^^jXfM, (3-24) 
where E a/3X ^ is the tensor of elastic moduli, i.e., 



tE 



2(1 + 1/) 



2v 



■a al3 a XfM 



(3.25) 



t = 2e is the shell thickness, and a a p = 5p. Setting 2e = t and using the relations (|3.21j) . we find 
that (|3.19|) and H3.25j) are equivalent. In matrix form, the right Cauchy-Green deformation tensor is 

C = F T F, 

where Vx = F is the deformation gradient and the Cauchy-Green strain (7) is 

G = ±(C-I). 

Assuming an isotropic film and the linear stress-strain relation in (|3,24|) . the Second Piola-Kirchoff 
stress can be written in matrix form as 



S = r 



where 



and the 2x2 cofactor matrix is 



G + ^Cof(G) rj 



tE 



(3.26) 



(3.27) 
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G and S are symmetric and by the spectral representation theorem, we have 

G = c^ni ® ni + <5 2 n 2 (g> n 2 , 
S = faxi! <g> ni + /x 2 n 2 ® n 2 , 

where ni and n 2 are orthonormal vectors. The eigenvalues of S (denoted by \x\ = t{5\ + z/J 2 ) and 
fj-2 = r(<5 2 + vSi)) are the principal stress resultants and the eigenvalues of G (denoted by 5i and 



5 2 ) are principal strains. Because we have assumed a linear stress-strain constitutive relation and an 
isotropic film, S and G have the same principle axes. The eigenvalues of C are the Cauchy strains 
which are denoted by A 2 where S{ = |(Af — 1) and <5j > —\. The him strain density is given by 

W f = iS : G. (3.28) 

In terms of the Cauchy-Green strains 5i, the standard membrane strain energy is given by 

W f {S 1 ,S 2 ) = (V + 5 2 2 + 2^82) , (3.29) 
2(1 v ) 

while in terms of the Cauchy strains, we have 

WfiXt, A 2 ) = |r (Ax 4 + A 2 4 + 2^A X 2 A 2 2 - 2(1 + v) + A 2 2 ) + 2(1 + ufj . (3.30) 

Next, we derive a few estimates that will be needed at a later time. From the definition of |F|, 
we have 

|F| 2 = Xj + Xl 

It follows that for < v < 1 

Ai 4 + A 2 4 + 2Ai 2 A 2 2 > A1 4 + A 2 4 + 2zvAi 2 A 2 2 > u{\\ + A|) + 2^Ai 2 A 2 2 = ^|F| 4 . (3.31) 
From (|3.3UI) and ()3.31() . we can hnd a lower bound for Wf, i.e., 

W f > |r(i/|F| 4 -2(l + ^)|F| 2 + 2(l + ^)) . (3.32) 
Applying l|3.31|) and ()3.32|) . we are led to 

|r (i/|F| 4 - 2(1 + ^)|F| 2 + 2(1 + 1/)) < W f < |r|F| 4 + |r(l + v). (3.33) 

We conclude this subsection with a lemma that will be useful in obtaining our existence results. 
Lemma 2 Let ki,k 2 > 0. 
(i) There exist a, p, 7 such that 

a\u\ 4 — p < Ki|n| 4 — k 2 |-u| 2 < -y I ii| 4 , (3.34) 
< a < n\ < 7, and p = k^/^ki — a). 



(ii) There exist constants 7, p such that 

Ki\u\ 4 + K2M 2 < l\ u \ A + P? (3.35) 
< ki < 7 and p = ^2/4(7 — Ki). 

Part (i). The second inequality in ()3.34j) is obvious. If — p = inf u {(«;i — a)|u| 4 — K2|"u| 2 } > —00, then 
we are done. Otherwise, it's easy to show that g(x) = (k\ — a)x 4 — K2X 2 has an absolute minimum 
at xq = \J ^2/2(^1 — a) and g(xo) = —K^j\{K\ — a). The first inequality in ()3.34[) follows with 
—p = g(xo). The proof of Part (ii) follows the proof of Part (i). 

3.2.3 Relaxation of the film strain energy density 

The energy density in (|3.28|) can lead to states where p\ or p2 are negative, corresponding to a 
compression. While this is reasonable for certain types of shells, the balloon film cannot support 
such a compression. Instead, the film will form folds or wrinkle. To tackle the problem of negative 
compressive stresses, we follow the methods introduced by A. C. Pipkin (see [24 ). In the following, 
let <5i and 82 be the principal (Cauchy-Green) strains for a typical facet T in a triangulation of 
Let pi and P2 be the corresponding principal stress resultants. In Pipkin's approach, a membrane 
M is decomposed into three distinct regions: 

S - slack region (5i < 0, 62 < 0), 

T - tense region (^1 > 0, p2 > 0), and 

U - wrinkled region (U = M\ SUT). 

We apply this classification scheme to each T in our triangulation of £1. On S the strain energy is 
assumed to be zero and on T the relaxed strain energy density is exactly the same as the standard 
strain energy density. On the region U, a modified Cauchy-Green strain G* is introduced. If G is 
the usual Cauchy-Green strain, then 

G* = G + /3 2 n®n, 

where n and t are (unknown) principal stress directions based on G*. Pipkin refers to — /3 2 n ® n 
as the wrinkling strain and G* as the elastic strain. The elastic strain is thought to represent the 
straining in an "averaged" wrinkled surface and leads to uniaxial stress on U in the form: 

S* = p t <g>t, 



where /i > and t is a unit vector orthogonal to n. For our exposition, we assume that t is the 
tensile direction, and n is a unit vector orthogonal to t. The parameter /3 2 and n are chosen in such 
a way that the following conditions are satisfied. 



n-S*n = 0, 
n-S*t = 0. 



For an isotropic material, S* can be written in the form: 



S + t/3 2 (n <g> n + z^Cof(n <g> n) T ) 



and it follows that 

p 2 = --n-Sn. (3.36) 
r 

If Wf = Wf(G), its relaxation is WJ = Wf(G*) where G* uses (3 2 from (|3.36|) and t is the principal 
direction that corresponds to a positive principal strain. One can show that on U, the principal strains 
of G* are in the form {62, —u5 2 } or {61, —v5\} (see ^E])- The principal stresses of S* = S(G*) are 
given by 



0, 5i < and 5 2 < 0, 

tE52^2 &> n2, /xi < and 62 > 0, 

tE5iiii (g) ni, [12 < and <5i > 0, 

/xini (g> ni + /i 2 n 2 (8) n 2 , /xi > and /i 2 > 0. 



(3.37) 



Wrinkling is modeled by replacing Wf with Wt where 



WK6 1 ,S 2 ;t,u,E) 



( 

l 

2 



Hes? 



2- 



\tE5\, 



tE 
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<5i < and <5 2 < 0, 
/ii < and #2 > 0, 
/i 2 < and 5i > 0, 

2-(5 2 + <5 2 + 2i^i<5 2 ), Mi > and /i 2 > 0. 



(3.38) 



In order to obtain existence results, we will need upper and lower bounds on WJ in terms of 
|Va;|. Since WJ is the largest convex function not exceeding Wf, we have 



WJ < W f . 

To obtain a lower bound for WJ, we consider the following. Let 

A = {(61,62) I <5i > -1,52 > -§} 
Ai = {(61,62) 6 A I 62 > -5 1 /v,5 2 > -v5i}, 
A 2 = A\Al 



(3.39) 
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Figure 5: (61,62) - domain for Wf and Wf 



See Figure [5J While Ai is unbounded, Wf = Wf for (81,82) € Ai. Since A2 is compact, there is a 
constant d such that 

d= max {Wf(8 1 ,8 2 )-WU8 1 ,8 2 )}. 

Wf(8±, 62) has a critical point when 5± + u5 2 = and 5 2 + v8\ = 0. Thus, (^1,^2) = (0, 0) is the only 
critical point. Hence, d is achieved on the boundary of A2 at a corner point. By inspection, one finds 
d is achieved at (-0.5, -0.5), WJ(-0.b, -0.5) = 0, and d = \r(l + v). Thus, 

Wf < Wj + \t(1 + v), (<5i,<5 2 ) G A. (3.40) 

Combining lj3~35|) - lpPn|) . yields 

W> - Jr(l + v) < Wf < Wf. (3.41) 

We conclude this subsection by noting if / is quasi-convex, so is / + C for any C £ K, for if 

/ / (u q ,vq, x ,£ + Vx(y)) dy > \D\f(u ,v , x ), 
Jd 

for every cube Dcfl, for every (u , v , x , ^ ) e!Jx M m x M nm and for every x G W 1,oo (D; M m ), 
then 

/ (/ («0, v , x , £ + Va:(y)) + C)dy> \D\ (f(u , v , x ) + C) . 

JD 

3.3 Constraints 

3.3.1 Tendon constraint 

Tendons are included in our balloon model through the use of constraints. We assume that a tendon 
is encased in a sleeve that is attached along the seam between two adjacent gores. We assume that 



the tendons are connected to a common point at the top of the balloon and to a bottom end-plate. 
The top of the balloon is free to move up and down the z axis and the bottom end-plate is fixed at 
the level z = 0. 

Remark 4 The interaction between the tendon and film is very complicated. Even without tacking 
to the sleeve, a tendon will not slide freely within its sleeve when a pumpkin balloon is fully inflated. 
The same can be said for a zero-pressure balloon which uses a load tape, a one-dimensional structural 
element that serves the same function as a load tendon. For simplicity, we have used the same tendon 
model for both the ZPNS and the pumpkin balloon. The assumption of constant tension enables 
us to include stiff tendons in the present model and analyze equilibrium configurations of a strained 
fully inflated shape. 

We first consider a single tendon in the deformed state x(T,i) for x £ T>. By the Sobolev 
Embedding Theorem (see, e.g., |20| Theorem 7.26, p. 171]; for a manifold setting, see PJ Theorem 
2.34, p. 55]), W 1,4 (n) is compactly embedded in C°>P(n) for any < (3 < 1/2. Hence, x £ X 
(or X g ) implies that x € C 0, ^(O,]R 3 ), and in particular, cti(v) = x\^ t is well-defined and Holder 
continuous. The length of sc(Ef) is denoted by L[£j]. A tendon is under constant tension, and the 
strain in the ith tendon is 

e i = (L\L i ]-L t )/L t . (3.42) 

In our numerical model, L[Sj] is approximated by 

n 

L[Ei] « \(*i(vj) - ati(vj-{)\ 

3=1 

where {vj} is a partition of [0, L c ], i.e., = vq < v\ < ■ ■ ■ v n -\ < v n = L c . 

The strain energy in a tendon of length Lt with stiffness K is St,% = St[T,i] = ^Kei 2 Lt. K has 
the units of force. Since the tendon cannot support compression, we utilize the relaxed tendon strain 
energy 

o, £i <o. {6A6) 

Since we assume that the tendons are inextensible, we set K = 1 and impose the conditions 

St fi <0, i = l,2,...,N, (3.44) 

which can be included as local constraints in the form of (|3.4|) . Because we are not using Va: to 
calculate the arc length of tendons, a continuity argument similar to the one used in Section |31 



shows X g is closed when g < includes a constraint in the form (|3.44|) . Consider (|3.7j) subject to 
constraints in the form of (|3.44|) . Let a; be a minimizer of this system where Aj, % = 1, . . . ,N are the 
corresponding Lagrange multipliers. We can interpret the Aj's in 

N 

E T (x,Vx)+Y,\Sl i . (3.45) 

i=l 

as the tendon stiffnesses that are needed to maintain equilibrium. In our numerical simulations 
presented here, we will consider cyclically symmetric solutions so only one tendon constraint is 
needed. Note, tendons are modeled in a different fashion in but our results here are consistent 

with those works when a very large value of K is used. We carried out numerical simulations when 
()3.44j) is replaced by an equality constraint and obtained essentially the same results as the '<' case 
with ei ~ 1CP 10 m/m. If we drop the requirement that the tendons are inextensible, we can replace 
the '<' in (|3.44j) with '>'. In this case, we find > for the solution where ei ~ 0.005 m/m. 

3.3.2 Volume constraint 

In an open system such as the ZPNS design, the balloon is open to the atmosphere at z = 0. This 
means that the volume of the balloon at equilibrium will adjust itself accordingly. In a closed system, 
the volume is fixed, and po is determined. Using the divergence theorem, the volume of the balloon 
S can be expressed as 

V(x)= [ 1 dV =\ [ xndS = [ f v (x,Vx)dA (3.46) 
Jd Js Jn 

where n dS = adj 2 Va; dA and fy(x,Vx) = gy(x,&dj 2 Vx), 

g v {x,A) = \x- A. (3.47) 

The same argument used to show fp is quasiconvex, shows fv( x , •) is quasiconvex. (|3.46j) can be 
evaluated exactly when S is a faceted surface. The volume constraint is in the form, 

V{x) -uj = 0, (3.48) 

where u>o is the target volume. 

In a closed system where (|3.48|l holds and po is unknown, we can assume —P{z) = bz. After 
finding an equilibrium solution by minimizing (|3.7|) subject to (|3.48j) . we find that the Lagrange 
multiplier A in 

£ T {x, Vx) + X(V(x) - uj ) 



yields the appropriate constant pressure term that is needed for equilibrium. Although we use 
Lagrange multipliers to handle constraints in our numerical approach, our theoretical model incor- 
porates such a constraint directly into the underlying solution space. 

Before we demonstrate that the volume constraint can be included in the definition of X s , we 
will need an estimate for |adj 2 A — adj 2 B|. 

Lemma 3 Let A, B G M nrn . There exists a constant (3 > independent of A,B such that 

|adj 2 A - adj 2 B| < 0(\A\ + |B|)|A - B|. (3.49) 

When nm = 6, $ = 3. 



We establish the proof for the case nm = 6, but it is straightforward to modify the argument for gen- 

(xi x 2 \ 
x 3 x 4 G M 3 ' 2 , consider X= [xi,x 2 ,...,a; 6 ] T . Define f(X) = [/i(X),/ 2 (X),/ 3 (X)] T 
x 5 x 6 J 

adj 2 X where /i(X) = x 3 x 6 - x 4 x 5 , / 2 (X) = -(xix 6 - x 2 x 5 ), / 3 (X) = x x x 4 - x 2 x 3 . For /i, we see 
Dfi(X) = [0, 0, X6, — X5, — X4, x 3 ] T . Since D 2 fi(~K) is independent of X, we write Z) 2 /i(X) = Hi, 
where 

/000 0\ 

1 
-1 
-1 
\ 1 / 

One easily obtains /i(X) = iX T HiX. For i = 2,3, we find D 2 /j(X) = Hi and /j(X) = iX T HjX, 



Hi 



where 



H, 
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Using the properties of f(A) and f (B), we obtain 
|f(A)-f(B)| 
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|(/i(A) - /i(B), / 2 (A) - / 2 (B), / 3 (A) - / 3 (B))| 
(A T Hi(A - B) + (A - B) T HiB, 
A T H 2 (A - B) + (A - B) T H 2 B, 
A T H 3 (A - B) + (A - B) T H 3 B) 



< \ ^|A T H i (A-B) + (A-B) T H i B| 2 ^J 

< \ ^|A T H i (A-B) + (A-B) T H i B|^) 



1/2 



= \Y^m (|A| + |B|)|A-B| 

1=1 

< /3(|A| + |B|)|A-B|. 
Since |Hj| = \/4 = 2, the proof is completed with (3 = 3. 

Lemma 4 Lei V(x) be defined as in Let \x k — sc 1 1,4 — ► where x,x k £ X and V(xk) = <^o> 

lim V(aifc) = V( lim a;*,) = V(a;) = ojo- 

fc— »oo k^oo 



For each A;, we have 



|F(sr)-wo| = \V(x)-V(x k )\ 

= U(x. adi2 Vx-x k . adh Vx k )dA. 
Jn 

Adding and subtracting x ■ adj 2 Va?fc to this last expression and applying Lemma El we are led to 



\V(x)-u \ 



= \ I (x — x k ) ■ adj 2 Vajfc + x ■ (adj 2 Va: — adj 2 V 'x k )dA 
Jn 

< / \x - x k \\Vx k \ 2 dA + 
Jn 

/ \x\(\Vx\ + \Vx k \)\Vx-Vx k \dA. 
Jn 



(3.50) 



Since x, x - x k ,\7x - Vx k G L 4 (0,M 3 ), Va?, Vx k G L 2 (ft,IR 3 ), and \Vx k \ 2 G L 4 / 3 (ft, IR 3 ), we can 
apply the generalized Holder inequality to (|3.5()|) and obtain 



\V(x)-u \ < 



\x — x k \ 4 dA 



\x\ 4 dA 



1/4 



1/4 



8 \ 3 / 4 

\Vx k \tdA] + 

1/2 



J \Vx\ 2 dAj + 



J \Vx k \ 2 dAj 1 (J \Vx - Va; fe | 4 dAj 



s 1/4 



(3.51) 



Since a^ — > a; G W l,i (p,, IR 3 ) and ||Va;fc||£4 is uniformly bounded, from 1)3. 5 lj) we see that \V{x) 
wo I = and the proof is complete. 



Remark 5 From Lemma 0J we see that V{x) — u>q = is a global constraint satisfying ()3.3[) in 
Lemma[I] When V(x) — ujq = is included among the constraints, g < 0, we denote the correspond- 
ing closed subspace of X by Xgf^ i n order to emphasize the dependence on ojq. 

4 Existence Results 

In the following, we write u = (u, v) £ £1. The total energy of the balloon system is taken to be the 
sum of the hydrostatic pressure potential and the film strain energy, 

£t(x,Vx) = J (Wf{u,x,Vx) + f P (x,Vx)^ dA. (4.1) 

Combining (|3.17|) and (|3.41j) . we find 

W f - \t{1 +u)- [\bR 2 + \p \R) |F| 2 < Wf + fp<W f + (±bR 2 + \p \R) |F| 2 . (4.2) 

Next, applying the appropriate upper and lower bounds on Wf established in (|3.33|) . we have 

\t (^|F| 4 - 2(1 + z,)|F| 2 + 2(1 + i/)) - \t{1 + v) - (±W? 2 + \ Po \r) |F| 2 
< W* f + f P < |r|F| 4 + ir(l + v) + (±W? 2 + \po\R) |F| 2 . 

Simplifying, we have 

Itv\F\ 4 - \ (r(l + v) + 2bR 2 + 4\po\R) |F| 2 

< W; + f P (4.3) 
< |r|F| 4 + ir(l + u) + {\bR 2 + \ Po \R)\F\ 2 . 

Let K\ = \tv and k 2 = j(t(1 + v) + 2bR 2 + 4\po\R) and apply Lemma Efi) to the first inequality 

in (|4.3[) . Choose ot\ so that < a\ < ^tv and define ui by the relation a\ = \tv\. We see that 

< v\ < v and by Lemma [2fi), we have 

|r^i|F| 4 - pi < |r^|F| 4 - ±(r(l + 2bR 2 + 4|p | J R)|F| 2 , 

where 

= {T(l + u) + 2bR 2 + 4\p \R) 2 
Pl ~ St{v - ui) 

Turning to the second inequality in l|4.3|) . and applying Part (ii) of Lemma El we choose 9\ > 1, and 
it follows that 

±r|F| 4 + (±W? 2 + |po|fl) |F| 2 < |r^|F| 4 + p 2 , 

where 

= (6^ 2 + 2|p |^) 2 
92 2t{6 1 - 1) ' 



Thus, 



and 



It follows that 



i^i|F| 4 - pi < W* f +f P < ir^|F| 4 + P2 + ir(l + u), 



|r^|F| 4 < W} + f P + pi< |r^|F| 4 + |r(l + i/) + p x + p 2 . 



(4.4) 



a;, Vx) = Wjf (w, a;, Vx) + / P (w, x, Vx) + pi (4.5) 
is quasiconvex and satisfies some additional properties that we summarize in the following lemma. 

Lemma 5 If f^(u, x, Vx) = Wj(u, x, Va?) + fp(u, x, Vx) + p\ and \x\ < R, then /y is quasiconvex 
and 

(i) 7 |A| 4 </^( M , a; ,A) <<5(1 + |A| 4 ) 

(ii) \ft(u,x,A)-f*(u,y,B)\ </?(l + |a;| 3 + |y| 3 + |A| 3 + 

|B| 3 ) (|x-y| + |A-B|). 



Property (i) is (|4.4I) with 5 = max{gT#i, |r(l + i/) + pi + p 2 } and 7 = |ti/i. We will show 
Property (ii) as follows. Consider 



f}(u,x,A) - ft(u,y,B) 



Wf(u,x,A)-Wf(u,y,B) 
+ (fp(u,x,A) -f P (u,y,B)) 



Since Lemma \A. II applies to Wf, there exists 0i > such that 

|W7(u,x, A) - W}(u,y,B)\ < ^(1 + |A| 3 + |B| 3 )|A - B|. 



(4.6) 



(4.7) 



To complete the proof of Property (ii), ±[|&(y • k) 2 + po(y ■ k)]k • adj 2 A is added to the last term 
in ((HD. Using ^J^-^J^j, we find 



f P (u, x, A) - f P (u, y, B)=I 1 + I 2 + I 3 + h, 



(4.8) 



where 



h 
h 
h 



(x • k) 2 - (y ■ k) 2 j k • adj 2 A, 
Po [{x • k) - (y ■ k)] k • adj 2 A, 
(y • k) 2 l k • (adj 2 A - adj 2 B) 



h = -po [y ■ k] k • (adj 2 A - adj 2 B) . 

By Young's inequality, we have 3|sc| |A| 2 < \x\ 3 + 2|A| 3 , and 3\y\ |A| 2 < \y\ 3 + 2|A| 3 . Using these 
estimates and Lemma El with (3 = 3, we find 

< |6|cc-y] (|a;| + |y|) |ad.j 2 A| 

< \b\x-y\ (\x\ + \y\)\A\ 2 

< ±6(j:r| 3 + |y| 3 +4|A| 3 ) \x - y\ 

< 2b(l + \x\ 3 + \y\ 3 + \A\ 3 + \B\ 3 )\x - y\. (4.9) 



Applying Lemma El with (3 = 3 and Young's inequality with 3|A| 2 < 1 + 2|A| 3 , we find 

1-^2 1 < bo | \x-y\ |adj 2 A| 

< 3bo| | » - y\ \M 2 

< bo| (l + 2|A| 3 ) \x-y\ 

< 2|po| (l + |a;| 3 + |y| 3 + |A| 3 + |B| 3 ) \x - y\. (4.10) 



Applying Lemma|3]with (3 = 3 and Young's inequality with 3|y| 2 |A| < 2|y| 3 + |A| 3 and 3|y| 2 |B| < 
2 1 2/ 1 3 + |B| 3 , we find 

|/ 3 | < i6|2/| 2 |adj 2 A-adj 2 B| 



< f%| 2 (|A| + |B|)|A-B| 



< \ (%| 3 + |A| 3 + |B| 3 ) |A-B| 

< 2&(l + |;c| 3 + |?/| 3 + |A| 3 + |B| 3 ) |A-B|. (4.11) 

Applying Lemma|3]with (3 = 3 and the estimates \y\ |A| < | + §|y| 3//2 | A| 3 / 2 < |(1 + \y\ 3 + |A| 3 ) 
leads to 

\h\ < bo I \y\ |adj 2 A- adj 2 B| 

< 3bo| \y\ (|A| + |B|)|A-B| 

< bo| (2 + 2|y| 3 + |A| 3 + |B| 3 ) |A-B| 
2bo| (l + |a;| 3 + |2/| 3 + |A| 3 + |B| 3 ) |A-B|. (4.12) 



< 



Combining (|4.9j l - (|4.12jl . we find 

\fp(u, x, A) - fp(u, y, B)| . , 

< 2(6 + bo|) (1 + W + \V? + I A| 3 + |B| 3 ) (\x - y\ + |A - B|) . 1 6 > 

Combining (|4.6|) . (|4.7|1 . and (|4.13|) . we see that Property (ii) follows with with (3 = 2max{/3i, 2(6 + 

bol)}. 



Since all of the forces are conservative, equilibrium is achieved at a minimum of the energy 
functional over the class of functions X g . The constant p\ serves only to increase WJ + fp by a 
constant, so the minimizer of Wj + fp is the same as the minimizer of WJ + fp + p\. We state two 
existence results for the equilibrium shape of an inflated balloon. 

Theorem 1 Closed balloon system. If \x\ < R, then 

OPg(^o)) inf|j(aj) = Jj^{u,x,Vx) dA \ x e X g(uJo) J 

admits at least one solution. 

Since we can always find an x £ xq + VFq 1 ' 4 ^; IR 3 ) parameterizing a spherical cap, we know that 
^g(wo) * s nonempty. In order to complete the proof, we need to show that f£ is quasiconvex and 
verify HI (i)-(iii) in Theorem lA.il Lemma shows that f T is quasiconvex and satisfies HI (i)-(ii) if 
| a; | < R. HI (iii) is satisfied because f T is independent of u and so r] = 0. Lemma 0] establishes that 
AgC^o) is closed. It follows that (Pg( wo )) has at least one solution. 



Remark 6 Eq. 1)1. lj) guarantees that the values of b and ujq are sufficient to lift the balloon system. 
If b = bd and loq = loq^ then the balloon envelope is fully deployed, and we are justified in ignoring 
the contribution of the film weight in the stress analysis. 

Remark 7 In general, one should not expect a unique solution of (Pg^)), especially if the balloon 
is not fully inflated. For example, if the lift generated by the gas bubble is not large enough to lift 
the balloon or the gas is compressed sufficiently (i.e., ujq < aj 0) d), then one expects to find multiple 
local equilibria with slack regions where a unique equilibrium configuration may not exist. For 
an ascending large scientific balloon, it is not uncommon to observe balloon shapes with a nearly 
periodic lobe pattern surrounding the gas bubble. These shapes are characterized by significant 
regions of excess folded material hanging beneath the gas bubble. Folded material can be handled 
by introducing constraints in the form covered by Lemma ^ and the theory developed in this paper 
applies to shapes with a period lobe pattern. 

Theorem 2 Open balloon system. If \x\ < R, then 

(P g ) inf j/O) = J ft(u,x,Vx) dA | x G X g | 

admits at least one solution. 



Description 


Variable 


Value 


Number of gores 


n g 


200 


Buoyancy (N/m 3 ) 


b 


0.068 


Tendon weight density (N/m) 
Film weight density (N/m 2 ) 


Wt 


0.094 


Wf 


0.344 


Payload (N) 


L 


4000 


End-plate diameter (m) 


di 


1.32 



Table 1: Shape finding parameters. For pumpkin balloons we assumed tb = 0.785 m and po — 200 Pa. For ZPNS 
balloons, po = Pa. Tendon weight and film weight are included in the shape finding process. 

Description Variable Value 

Film Young's modulus (MPa) E} 404.2 

Film Poisson ratio Vf 0.825 

Film thickness (/jm) t 32 

Table 2: Mechanical properties of polyethylene film. 

Since the subspace X s is closed by Lemma ^ the proof of Theorem |2] is the same as the proof of 
Theorem n with the space Xg^ ) replaced by X s . 

5 Numerical results on strained balloon shapes 

In this section, we present numerical solutions of our model. We consider a zero-pressure natural 
shape design and a pumpkin design. We assume the balloon film is 32/mi polyethylene and the 
balloons have no caps. We use the same tendon weight density and a suspended payload of 4000 N 
for both balloons in shape finding. The specific buoyancy at float is b = 0.068 N/m 3 and corresponds 
to an altitude of roughly 35 km. For the pumpkin balloon, vb = 0.785 m and po = 200 Pa. See 
Table ^ for a summary of shape finding parameters that were used. We used mechanical properties 
that were determined by the Balloon Lab at NASA's Wallops Flight Facility for 0.8 mil polyethylene 
film at room temperature (this is comparable to the temperature at float during nominal daylight 
conditions). For an isotropic three-dimensional material, one can show that Poisson's ratio is between 
and 0.5 (see, |13l p. 129]). However, in lab experiments using cylinder tests, one finds that for 
a thin polyethylene film, Poisson's ratio is greater than 0.5 (see, e.g., ^^)- We extrapolated the 
findings for 0.8 mil to 1.5 mil film to estimate Young's modulus for our simulations. Mechanical 
properties are summarized in Table El 

Next, we summarize our findings when we numerically solved Problems (P g ) and (P g ^ ))- In all 
cases, we assumed a cyclically symmetric shape and so we needed to solve for one-half a gore. We 
divided a half-gore into three strips with 100 triangles per strip. We used our own finite element 



code written in Matlab to compute the energies and constraints as described in Sections 3-5. We then 
used fmincon from Matlab's Optimization ToolBox to solve a constrained minimization problem where 
was the objective function and g < were the constraints. Gradients for /y,g, and the Hessian 
of /j. were computed analytically. When a nonlinear constraint (i.e., volume or tendon constraint) 
is imposed, we used a "medium-scale" sequential quadratic programming quasi-Newton's method; 
when there were only linear constraints, a "large-scale" projected trust-region Newton's method was 
used (see, |14| for further details on these algorithms). 

We discretized the three-dimensional surface produced by the shape finding process and used 
this as an initial guess to start the solution process. For each of the scenarios considered, we present 
plots of the "averaged" principal strains and "averaged" principal stress resultants for each adjacent 
pair of rectangles in a strip. If a meridional strip has 2M triangles (numbered from top to bottom) 
and fi2i-.i i and ^2i,i are the meridional principle stress resultants, then we plot \{fi2i~i,i + M2i,i) 
for % = 1, . . . M. Note, the directions for fi2i-l,l an d M2i,l need not coincide, but this measure seems 
to work satisfactorily for summarizing the data. A similar convention is followed for averaging the 
circumferential stress resultants and the principal strains. 

5.1 Zero-pressure natural shape balloon. 

We considered two scenarios involving a ZPNS balloon. 

(a) Closed system V{x) = ujq. In this case, we considered a closed system with a target volume 
u>o = 137,023 m 3 . The averaged principal strains are plotted in Figure Efa) and the averaged 
principal stress resultants are plotted in Figure Efb). The 81 and 82 directions corresponds 
roughly to meridional and circumferential, respectively. We found that the maximum principal 
strains were less than 0.03%. The meridional stress resultants are plotted in the top graph 
in Figure EJb) and increase monotonically from the bottom of the gore to its top. From 
Figure EJb), we see that near the top one-third, the film is in a biaxial state and in the middle 
one-third it is in a uniaxial or near uniaxial state (i.e., fi2,i ~ 0). 

(b) Open system po = 0. In Section IS~lT b). we considered the same conditions as those in Sec- 
tion E^a), except we dropped the volume constraint and set po = 0. The results were nearly 
identical to those obtained in Section 15.1( a) and for this reason plots were not included. We 
found that the difference between the volumes in Section fo.K a) and (b) was less than 10 -5 m 3 . 



In the following, principal strains and stress resultants are to be interpreted as 'averaged' principal 
strains and 'averaged' principal stress resultants, respectively. 

5.2 Pumpkin balloon without tendons. 

To get a sense of the important role that tendons play in a pumpkin balloon, we computed a solution 
for a pumpkin balloon without tendons. Since it was clear that a single layer of 32 /im film was not 
sufficiently strong, we quadrupled the thickness for this numerical experiment. From Figure Ufa) we 
see that the strains are nearly 7.0% and the principal stress resultants are quite large (a maximum 
of 8 kN/m). It is interesting to note that near the ends, the deformed gore is in a state of biaxial 
tension, where in the middle of the gore the tension is uniaxial and parallel to the gore length. In 
reality, the film is a viscoelastic material and under the right conditions may be able to maintain the 
integrity of the gas barrier at such a high strain. However, the peaks in Figure [7] suggest this loading 
scenario should be avoided. 

5.3 Pumpkin balloon with tendons. 

In this case, we considered a single layer of 32 ^im polyethylene film and added inextensible tendons 
via a constraint (i.e., (|3.44j0 . In Figures Eta)-(b), we present the principal strains and principal stress 
resultants. We see that when tendons are added, the maximum principal strains are reduced from 
7.0% in Section [5.21 to 1.5%. The maximum principal stress resultants are reduced from 8 kN/m in 
Section f5. 21 to 225 N/m. It is interesting to note from the plot of the circumferential stress resultants 
(bottom graph of Figure Etb)), the film is in a state of biaxial stress along the center of the gore. 
But near the equator and the sub-tropic regions, as one gets closer to the tendon, the film is in a 
uniaxial state. This makes sense because one would expect there to be wrinkling in these locations. 

5.4 Pumpkin balloon with shortened tendons 

In Figures IHKa)-(b), we plot the principal strains and principal stress resultants when shortened 
tendons are utilized. Basically, we reduced the tendon length in Section 15.31 by 2%. Overall, we 
see that the principal strains and principal stress resultants are reduced. However, the maximum 
principal strains and principal stress resultants are not reduced by a significant amount. Comparing 
Figure IBth) and Figure EJb), we see more wrinkling in Figure Efb), especially in the area adjacent 
to the tendons. 



While there are many other factors that go into the design of a gore, a balloon of the type discussed 
in Section 15.41 is not very efficient. Ideally, one would like to have the gore uniformly loaded. The 
presence of wrinkles in the fully inflated shape suggests that there is excess material that is doing 
little to carry its share of the load. Furthermore, it has been demonstrated that too much tendon 
shortening can be detrimental, leading to instability in the strained cyclically symmetric float shape 
and impeding proper deployment (see, [H]-[Z|)- 

6 Conclusions 

Motivated by the problem of a large scientific balloon at float altitude, we presented a mathematical 
model for a strained inflated wrinkled membrane loaded by hydrostatic pressure and constrained 
by tendons. We assumed the balloon is constructed of a thin linearly elastic isotropic material. 
Load tendons are included through the use of constraints. Using direct methods in the calculus 
of variations we established rigorous existence theorems under general conditions. We computed 
numerical solutions based on our model and estimated the principal strains and principal stress 
resultants under nominal loading conditions at float altitude. Our theoretical results establish a solid 
foundation for our mathematical model and affirms the use of numerical computations to estimate 
film stress resultants, information that is valuable to the balloon designer. An efficient balloon design 
keeps weight to a minimum, but the balloon must be of sufficient strength to operate safely over 
its service life. Analytical and computational tools are demonstrated here that can help the balloon 
designer balance the competing factors of a stronger, albeit heavier balloon versus a lighter more 
efficient one. 
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A Direct methods in the calculus of variations 



For the convenience of the reader, we record two results by Dacorogna. 

Lemma A.l Dacorogna ]lb\ Lemma 2.2, p. 156] Let f : IR" — > IR be convex in each variable and 
let 

|/(x)|<a(l + |af) 
for every x £ M n and where a > 0,p > 1. Then there exists (3 > such that 

\f(x) - f(y)\ <p(l + Ixr 1 + lyr 1 ) |x - y\ 

for every x, y G IR n . 

Theorem A.l Dacorogna \16\ Theorem 2.9, p. 180] Let Vt C M n be a bounded open set. Let 
f : O X M m x M nm — > M be continuous and quasiconvex and satisfying 

(i) 7 | A| p < f{u, x, A) < 5(1 + \x\ p + \A\ P ); 

(ii) \f(u,x,A)-f(u,y,B)\<p(l + \x\P- 1 + \y\P- 1 + 

| A |p-i + |B|p-1) (|aj-y| + |A-B|); 

(iii) \f(u, x, A) — f(v, x, A) | < n(\u — v\) (1 + |a;| p + | A| p ) where n is a continuous increasing func- 
tion with n(0) = and p > l,a,(3, 7 > are constants. 

Let 

(P) inf|/(aj) = J f(u,x(u),Vx(u))dA : a; G x + W 1,p (O; iR m )| (1.1) 
then (P) admits at least one solution. 
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Figure 6: Zero-pressure natural shape balloon with ten- Figure 7: Pumpkin balloon without tendon constraints 
don and volume constraints. and nominal thickness quadrupled. 
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